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A model of the open Universe described by a conformally flat 4-metric is considered. The gravitational equations 
with a perfect Pascal fluid as a source are reduced to the nonlinear equation of oscillations. It is proposed to consider 
this equation as the Mathieu equation encompassing some cosmological models. It is shown that the behaviour of 
the Universe state function is in qualitative accordance with the Big Bang scenario. 
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The interest in cosmological models does not fall in 
the recent decades, in particular, due to multiple obser- 
vational data, which in accordance with perfecting the 
instruments and methods of observation are constantly 
improved. According to the accepted modern views, the 
matter filling our Universe changed its state during the 
evolution: from a stage of physical vacuum to an ultra- 
relativistic stage and further to the Friedmann stage. It 
means that topical in the theory is a search for exact 
solutions of the gravitational equations with the matter 
equation of state depending on a space-time point. 

Let us consider a model of the open Universe of 
Friedmann type, which can be described by the con- 
formally flat metric 
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- r is a square 
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The Einstein gravitational equations are written as 



a/3, 



(2) 



with the Einstein tensor G a p and the energy-momen- 
tum tensor (EMT) of the Pascal perfect fluid 



T a fs = e ■ u a up + p-b a 



/9> 



(3) 



where e is the energy density, p is the pressure, b a p — 
(u a U/3 — g a p) is a projector onto the 3-space, g a p is 
the metric 4-tensor, u a — exp(cr) ■ S :Ct is the 4-velosity, 
Greek indices take the values 0,1,2,3. The velocity 
of light and the Newtonian gravitational constant are 
chosen to be equal to unity. 

Substituting the metric (1) into Eqs.(2), we obtain 
two differential nonlinear equations (for the energy den- 
sity and pressure) 
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dy dy 
dx ax 
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(4.1) 
(4.2) 



d 2 y _ _]f 
dx 2 ~ 4x 4 ' 

where the new function y = e a / 2 and the new variable 
x = 1/S have been introduced. 

The last equation can be reduced to a representation 
of Newton's second law 



dx 2 



= F 



(5) 



for the potential force F = —dll/dy and the potential 
function U = £l 2 (x)y 2 /2, if we define the right-hand 
side of Fq.(4.2) as the force F [1]. 

Thus we have the differential equation of nonlinear 
oscillations 



J + n>(x). y = o, 



(6) 



where the function fl 2 (x) is the squared frequency. 

To solve the problem of finding a solution to this 
equation it is suggested to reduce this equation to the 
Mathieu equation, extending some special cases ob- 
tained earlier, 

^l + B 2 [l + hcos{ 1 x)]-y = Q, (7) 

where B, h, 7 are constants. 

If ft 2 — 0, i.e. the pressure is absent (p = 0), 
we obtain the Friedmann solution for the open Universe 
filled with incoherent dust, 



y F = 1 - A3 



(8) 



where A = const is determined the by modern value of 
dust density in the Universe 
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(9) 
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When fl = B 2 = const , we have an axact cosmolog- 
ical solution [2] for the open Universe filled with dust 
and electromagnetic equilibrium radiation, 



y BS = ^l + {A/B fcos{Bx - ao), 



(10) 



where tgao = A/B, the constant value B is connected 
with the density of equilibrium cosmological electromag- 
netic radiation 



B 2 
12^, 



(11) 
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and the solution (10) at large times (S — > oo; x 
passes through the Friedmann solution (8). 

The Mathieu equation (7) may be rewritten in a 
dimensionless form by introducing the variable C = 
B/S = Bx : 



d 2 y/d( 2 + [1 + hcos(j(/B)} -y = 0. 



(12) 



For solutions of the Mathieu equation ther are reso- 
nance regions in accordance with the values of the pa- 
rameter 7 = 2B/n, n = 1,2,... Here we shall choose 
the first resonance region (7 = 2B). Then Eq.(12) will 
be rewritten as 



d 2 y/d( 2 + [1 + hcos(2Q] -y = 0. 



(13) 



In this case, according to Floquet's theorem the so- 
lution of the Mathieu equation may be substituted as 



y = e^ce 1 ((,h)-C-e-^se 1 ((,h), 



(14) 



where fi is a chracteristic parameter, cei and sei are 
Mathicu's functions having the form 



cei(C) 



k 



A k cos(2k + 1)C; 



sei(C) = ^A k sin{2k + l)C,. 



(15.1) 
(15.2) 



The Mathieu functions, the signs and constant C = 
fi + A/B are chosen here so that in the transition to the 
modern epoch (S — > 00; i^O) the solution (14) will 
pass through the Friedmann solution (8), and at /1 = 
and k = 1 it will be the same as the solution (10). 

By means of parameter selection it is possible to con- 
struct the physical behaviour of the function of state 

M x \ = P= 1 x ■ V ' d 2 y/dx 2 

e 3 (dy/dx) ■ (xdy/dx — y) ' 

which is at each instant S an equation of state of the 
Universe. 

It is clear that the first approximation of the solution 
(14) will be written as 

y w e MC cos ( - C ■ sin ( 
= e^ B ' s cos{B/S) - C ■ e-^ B ' s sm(B/S). (17) 



The numerical value of the parameter u = 0.56 
is found for Eq.(13) using a computational method of 
a characteristic exponent from Hill's determinant [3]. 
Other parameters may be selected so that we had the 
qualitative picture of model behaviour of the Universe's 
function of state according to the Big Bang scenario. At 
the first instants (when we have a cosmological singular- 
ity) the Universe, being initially in the state of physical 
vacuum (S = 0,(3 w — 1), in an inflationary phase is 
warmed up and reaches the radiation phase (ft w +1/3) 
and then, being expanded and cooled, slowly evolves to 
the modern epoch (/? w 0). 
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